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INTRODUCTION 


Stokes’  equation  gives  the  mathematical  relationship  between  the  geometric  shape  of  the  geoid  and  the  mag- 
nitude of  gravity  on  the  geoid.  The  equation  is 

N ’ <"  I 

where 

R = mean  radius  of  the  earth 

G = mean  value  of  gravity  over  the  earth 

S(i //)  = Stokes’  kernel 

The  integration,  extended  over  the  unit  sphere  using  a global  set  of  gravity  anomalies  Ag,  produces  the  geoid- 
ellipsoid  separation  at  a given  point. 

In  practice,  the  integration  is  not  extended  over  the  entire  globe,  but  only  up  to  a spherical  distance  i pQ  from 
the  computation  point  due  to  a lack  of  sufficient  surface  gravity  data.  The  error  5N  that  results  from  neglecting  data 
beyond  \pQ  has  been  analyzed  by  Molodenskii  (Reference  1)  and  is  described  in  Moritz  (Reference  2).  The  mag- 
nitude of  this  error  may  be  reduced  if  the  gravitational  potential  is  known  in  terms  of  a finite  spherical  harmonic 
expansion  derived  from  satellite  observations  by  applying,  for  instance,  the  computational  technique  of  Rummel 
and  Rapp  (Reference  3).  Equivalently,  the  error  may  be  reduced  by  removing  known  lower-degree  harmonies  from 
the  gravity  anomaly  data  used  when  Equation  (1)  is  applied  to  a spherical  cap  assuming  the  lower-order  harmonic 
components  of  N are  known.  In  addition,  one  might  also  consider  the  possibility  of  modifying  or  replacing  Stokes’ 
kernel  S(\p)  with  a function  which  produces  a better  estimate  of  N using  a limited  data  set.  The  purpose  of  this 
report  is  to  determine  what  decrease  in  5N  may  be  expected  if  harmonics  are  removed  from  S(i//)  or  Ag  individually, 
or  from  both  simultaneously. 


INFLUENCE  OF  DISTANT  ZONES  ON  STOKES’  EQUATION 


CASE  I:  STANDARD  CASE 

The  error  in  the  geoid  undulation  due  to  neglecting  gravity  anonaly  data  outside  a spherical  cap  a of  radius 
is  given  by  Moritz  (Reference  2)  as 


S(i p)  Ag  sin  d^  da 


(2) 


where  the  subscript  I refers  to  the  removal  of  harmonics  from  the  data  only  (Table  1).  Defining  a new  function 
S(tf)as 


1 


S(*)  = 


f 0 
S(^) 


0 < tit  < 


\pn  < \l>  < It 


Equation  (2)  becomes 


/n  fin 

. ,1,0 


da 


The  function  S(0)  has  the  following  spherical  harmonic  expansion: 


S(*)  = £ Qn  p„  (cos  *) 

n=0  1 n 


where 


3"  ~f  ‘ 

•'ft 


Qn  = / s(^)  pn  (cos  sin  \p  d\p 


■X 


S(0)  Pn(cos  sin  i//  di// 


Following  Moritz  (Reference  2),  the  RMS  error  due  to  neglecting  distant  zones  is  the  square  root  of 

o.2'» 

4G2  n=2 

where  it  is  assumed  the  zero  and  first-order  harmonics  of  the  gravity  anomaly  field  are  zero. 

Ag(0,A)  = £ Agn(0,A) 
n=2 

The  quantity  cr  is  the  anomaly  degree  variance  for  degree  n defined  by 

Cn  = E{A8n  } 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


where  E is  the  expected  value  operator.  The  anomaly  degree  variance  is  a function  of  the  spherical  harmonic  coeffi- 
cients of  degree  n 


II 

n = £ + _BL> 


m=0 


(10) 


where  the  bar  indicates  normalized  coefficients.  The  coefficients  Qn  are  known  as  the  Molodenskii  coefficients, 
and  much  is  known  of  their  properties. 
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CASE  II  REMOVAL  OF  M LOWER-DEGREE  HARMONICS  FROM  S(tf/) 


Equation  (1)  may  be  written  in  the  following  form: 


S(vl')  Ag  sin  t|/  di^'  da 


jl  r r 

4rrG  J , . , J _n 

l//  = l//  %/Q-  0 


S(i//)  Ag  sin  iM  ^ da 


The  kernel  S(i^)  has  the  spherical  harmonic  expansion 


Now  define 


S(^)  = E Pn(cos  *) 

n=  2 n 1 


sw>  = sw  - E Pn(c°s 


S(0)  = S(*)  + SU) 


Substituting  Equation  (13b)  into  Equation  (1 1)  gives 


N = + I2  + l3 


where 


I >-*- 
1 4rrG 


P P 

=0  ^ Q=0 


S(i^)  Ag  do 


£ £ Swae 


da 


£ £ s™,*,4Eio 


01) 


(12) 


(13a) 


(13b) 


(14a) 


(14b) 


(14c) 


(»4d) 
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Considering  the  definition  of  S|n(^)  and  Equation  (8),  the  integral  I3  becomes 

I3  = N2  + N3  + ...  + Nm  (15) 

which  are  the  first  m harmonic  components  of  N at  the  computation  point.  Assuming  I3  to  be  known,  the  error  due 
to  integration  over  a limited  cap  becomes 


L *** 


do 


Define  the  function  S as 


s(tio  = 


0<\p<\p 


0 


S(*) 


0 < ip  < \p 


0 


S<I^)  — Sm(^/)  \pQ^\p<n 


The  function  S has  the  Legendre  series  expansion 


where 


■/' 

*/0 

=/' 

0 

= Qn  - 


SOfr)  = £ Q Pn  (cos  *) 

n = 0 n n 


S(<p)  P„(cos  ip)  sin  \p  di p 

f 9 

S(^)  Pn(cos  ip)  sin  ip  d\p  — / S (\p)P  (coup)  sin  \p  d\p 

Jit 

0 

f" 

I Sm(^)  pn(cos  \p)  sin  ip  dtp 

J<b 


Therefore, 


SN„  - 


JL  f O' 

4”°A  X 

* f ['■ 

4ffG  Jo  Jq 


S(^)  Ag  do 


S(*)  Ag  do 


(16) 


(17) 


(18) 


(19) 
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--Sg  £<2"+i>5 


n=0 


/'/ 


2tt 


n Ag  Pn(cos  i/z)  sin  i//  d^/  da 

70  •'() 


R ~ 4ff  Ag„ 

8*G  2 (2n  + l)Qn  (2n+  1) 


2G 


E 0„ 


n=0 


The  RMS  error  is,  following  Moritz  (Reference  2), 


8nm  * j?  E «. 


CASE  III:  REMOVAL  OF  M LOWER-DEGREE  HARMONICS  FROM  Ag 


From  Equation  (8), 


Ag  = J2  Agn(0,X) 
n=2 


where  the  n'th  degree  hamonic  Ag  is  given  by 


A<?n  = E + 

m=0 


Removing  the  first  m harmonics  from  the  data  set  gives 


Ag  = Ag  - £ Agn 
n=2 


Equation  ( I ) becomes 


n-^/X 


2n 


S(i^)  [Ag  + y;  Agn]sin  \j/  dt//  da 
n=2 


= ll  + l2 


(20) 


(21) 


(22) 


(23) 


(24a) 
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iff 


S(  4/ ) Ag  sin  i//  d i//  da 


= 4^G  22  f f W)**n  sin*d*1 

n=2  •'0  •/() 


The  second  integral  reduces  to 


l2  = N2  + ...  + Nm 


which  is  assumed  known  since  Ag2 Agmare  assumed  given.  Therefore, 


N = N2  + . . . + 


Nm  * S5  f'f 

*'o  Jo 


S(\p)  Ag  da 


S(i^)  Agda 


Ttie  error  due  to  distant  zones  is  the  second  integral  in  Equation  (26). 


"I  47tCi 


iff 

*/0  */0 


S(^)Ag  da 


= l S(^)Agda 

= (2n  + ')Qn  / / AgPn(cos  ^)sin  t// d^  I 

n=0  ^0  *'o 


For  n = 0,  I , . . . , nt 


/TV 

*/0  *'o 


(cos  sin  \p  d^/  da  = 0 


since  Ag  contains  no  harmonics  of  these  degrees.  Therefore 


It  should  be  noted  that  this  error  is  the  same  as  that  obtained  using  the  Rummel  and  Rapp  (Reference  3) 
procedure  with  the  first  m harmonics  applied  (neglecting  here  corrections  for  atmospheric  effects  and  zero-order 
undulation  correction). 


CASE  IV:  REMOVAL  OF  M LOWER-DEGREE  HARMONICS  FROM  BOTH  S(\p)  AND  Ag 


Using  equations  (13b)  and  (23). 


Equation  ( 1 ) becomes 


where 


S(i//)  Agd  a 


S(t//)  Agn  do 


4 4rrG 


h £ IT 

n--  *'()  •'0 


Sn,(^)Agn  da 


(31e) 


By  orthogonality,  l7  and  Ij  vanish.  Also,  by  definition  of  S (Equation  13a |, 


4 47rG 


h£  £ fH 


do 


If  n ^ j,  the  partial  sum  is  zero  by  orthogonality.  When  n = j,  the  nonvanishing  terms  sum  to 

I = R 2n  + 1 47t 

4 4rrG  n - 1 2n  + I 


n = 2 


since 


f f P„(cos  i//)Ag  do  = -ff  ^Bn 

Jo  J0  2n  + 1 


Therefore, 


14  = N7  + . . . + Nm 


Thus,  Equation  (31a)  becomes 


N = N,  + . . . + N + ~~ 

2 4ttG  , 

•Ja  = 0 


271  ~ ~ 
S(ip)  Ag  do 


-S_  f'  f } 

4nG  / / 

•/<j=0 


2rr 

S(l//)  Agdo 


(32) 


(33) 


(34) 


Assuming  the  first  m harmonics  of  N are  known,  the  error  in  N due  to  neglecting  distant  zones  is 


6N 


IV  4trG 


k /'/■' 
•'*  */0 


S(0) Agdo 


8 


n 


.i  r f1' 

4*Gii  Jn 


S((//)  Ag  da 


n=0 


,/■'  f2’ 

J\h=0  •“'a=0 


Ag  Pn(cos  i/<)  do 


For  n = 0,  1 , . . . , rn  the  integral  in  Equation  (35 ) vanishes  giving 


6N,v  = £ £ QnAgn 
n=m  + l 


The  RMS  error  is  the  square  root  of 


«N?V  ~ £ ol  c„ 

4G  n=m  + l 


Tables  1 and  2 present  summaries  of  Cases  I through  IV  and  of  RMS  error  formulae,  respectively. 


Table  1 . Summary  of  Cases  1 through  IV 


(35) 


(36) 


(37) 


Case 


Equation  for  Undulation 


N ~~£cffs(*)Ag'‘°  * SN| 


"■".•-".’St// 


S(i/0  Ag  do  + 6Nn 


III 


N = N2  + ...  + Nm+^  / / S(*)Ag 


do  + 6N, 


til 


IV 


N = N2  + .. 


*Nn.-,4ic//  a*»5d«*5Nlv 


•-“V 


Table  2.  Summary  of  RMS  Error  Formulae 


r 

! 
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Case 

M Harmonics  Removed 

Square  of  RMS  Error 

I 

- 

D 2 oo 

— V02c 

4G2  “2  " " 

11 

S(*) 

R2  r’  pr2 

4G2  A-;0"  C" 
n=2 

III 

Ag 

— V Q2c 

4g2„4v  " 

IV 

S(H  Ag 

R2  v n2 

C" 

™ n=m  + l 

NUMERICAL  EVALUATION  OF  ERROR  FORMULAE 


The  RMS  error  for  each  case  was  computed  for  m equal  to  6, 1 2,  and  1 8.  In  these  computations  the  anomaly 
degree  variances  were  computed  using  Kaula’s  rule 


cn  = mrrj  mgals2  n = 3>4-  --  <38) 

and  c2  was  taken  as  10  mgal2.  The  constant  R2/G2  was  taken  as  42.3  m2/mgal2. 

The  Molodenskii  (Reference  1)  coefficients  Qn  were  computed  using  two  algorithms.  Initially,  the  Qn  were 
computed  using  an  algorithm  of  Hagiwara  (Reference  4);  however,  due  to  numerical  difficulties,  only  the  first  50 
terms  in  the  error  expressions  could  be  summed.  The  computations  were  repeated  using  the  algorithm  of  Paul 
(Reference  5).  An  increased  number  of  terms  (250)  were  summed  in  this  case. 

The  coefficients  Qn  were  computed  using  the  following  procedure: 

Qn  = Qn  - / Sm  Pn  (cos  1 1>)  sin  $ dj/  (39) 

0 


ij 


! 


; 


i=2  1 


(40) 


10 


- ->i  + 1 f ” 

Qn  = \ I Pj(cos  ^)Pn(cos  0 ) sin  4/  d\l/ 

j=2  A 


where  Qn  = Qn(^0) 

Qn  = Qn(*0> 


Taking  x = cos  \p 


P|  (cos  4/)  Pn  (cos  ip ) sin  \p  d\l/ 


/cos^ 

°Pj(x)Pn(x)dx, 


and  using  the  identities 


Pn(x)dx  = 2^n  (2n_,) 


P‘.,(X)  dx 


+ * fPn(X)  + Pn-.(X)I  - 2Pn-l  Pn 


/*b  (x(n -k)P  (x)P.(x)  + kP  (x)P.  ,(x)  - nP  ,(x)P.(x)] 

J Pn(x)Pk(x)dx  = 4 in-kXn+k+li  ' 


liquation  (41 ) may  be  evaluated. 


CONCLUSIONS 


Figures  1 through  4 give  the  root  mean  square  error  for  Cases  I through  IV  where  250  terms  have  been 
summed  in  the  error  equations.  In  Cases  II  through  IV,  the  first  m harmonic  components  of  the  geoid  undulation  N 
were  assumed  known.  An  examination  of  these  figures  shows  that  for  spherical  caps  of  radii  < 60  , removal  of 
known  lower-degree  harmonics  from  the  cap  data  alone  produces  the  least-expected  error  due  to  neglected  /ones. 
Removal  of  lower-degree  harmonic  components  from  the  kernel  function  S(i^)  does  not  minimize  the  expected  error 
for  cap  sizes  normally  considered.  Thus,  the  former  procedure  (Case  III)  is  recommended  it  the  assumptions,  upon 
which  these  error  formulae  are  based,  are  satisfied. 

Finally,  a comparison  of  the  results  based  on  50  and  250  terms  in  the  summation  demonstrated  that  satisfac- 
tory convergence  of  the  series  had  been  obtained. 


RMS 

ERROR  ( m ) 


Figure  4.  Root  Mean  Square  Error, 
Cases  II  through  IV,  m = 18 
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